Superfluid-spiral state of quantum ferrimagnets in magnetic field 
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We study the phase diagram of one-dimensional quantum ferrimagnets by using a numerical 
exact diagonalization of a finite size system along with a field-theoretical non-linear a model of 
the quantum ferrimagnets at zero temperature and its effective description in the presence of the 
external magnetic field in terms of the quantum XY-model. The low- and the high-field phases 
correspond respectively to the classical Neel and the fully polarized ferromagnetic states where in 
the intermediate magnetic field (/i cl < h < h C 2), it is an XXZ+h model with easy plane anisotropy, 
which possess the spiral (superfluid) states that carry the dissipationless spin-supercurrent. We 
derive the critical exponents, and then will study the stability of the XY spiral state against these 
spin-supercurrents and the hard axis fluctuations. We will show a first order phase transition 
from the easy plane spiral state to a saturated ferromagnetic state occurs at h — h C 2 if the spin- 
supercurrent reaches to its critical value. 

PACS number: 76.50.+g, 75.50.Gg, 75.10.Jm 



> 
OS 



o 
o 



i 

C 

o 
o 



I. INTRODUCTION 

Recently, antiferromagnetically coupled mixed-spin 
chains with an alternating array of two kinds of spins 
have attracted interest among researchers Jl|-^|. In- 
tegrable models of mixed-spin antiferromagnetic chains 
were constructed by de Vega and Woynarovich H and 
the simplest case of such chains with spins S = 1 and 
1/2 were subsequently studied Since these intc- 

grable models are exactly solvable, they are very useful 
for studying (quantum) statistical mechanical properties. 
Although ferrimagnetic spin chains exhibit both ferro- 
magnetic and antiferromagnetic features, they show some 
peculiar, and sometimes surprising, features uncharacter- 
istic of either the ferromagnet or the antifcrromagnet — an 
example being the existence of gapless excitations with 
very small correlation length. It is important to under- 
stand these features more clearly. In this paper we first 
present a numerical study on the quantum ferrimagnetic 
spin chain in magnetic field by an exact diagonalization 
of a finite size system. We show that the low- and the 
high-field phases correspond respectively to the classical 
Neel and the fully polarized ferromagnetic states where 
in the intermediate magnetic field (h c i < h < h C 2), it is 
an XXZ+h model with easy plane anisotropy an d 
the critical exponents are derived numerically. Then the 
physical properties of the quantum ferrimagnets in the 
presence of the uniform magnetic field, as described by 
a phenomcnological field theory based on a continuum 
non-linear a model (NLerM) formulation are examined 
further and will explore its novel features. When h > h&, 
the quantum ferrimagnets can be described by an easy- 
plane anisotropy state where the broken U(l) symmetry 
of spins corresponds to a topological spiral (superfluid) 
state. A mctastable state with a non-zero topological 



spin-supercurrent J s , is the superfluid state of the or- 
der parameter in XY-plane. Such possibility, may be led 
to the remarkable spin-dependent transport phenomenon 
in the ferrimagnetic chains, i.e., the spin-supercurrent is 
carried collectively (rather than by quasiparticles). Be- 
cause the spin current is non-zero when the system is 
in equilibrium, it flows without dissipation. The spiral 
state can be characterized by a wave- vector Q, i.e., the 
soliton (saw tooth) lattice spacing a — 2ir/Q is a length 
scale through which the in-plane phase (p of the order 
parameter changes by 2n. Hence, the spin-supercurrent 
is characterized by these wave-vectors, e.g., J s oc Q. In 
analogy with the dissipation mechanism of the supercur- 
rent in a one-dimensional superconductors ||, at large 
J s , hard-axis fluctuations of the order parameter becomes 
unstable and a first order phase transition to a uniform 
in-plane state (ip = cte) takes place, and the long range 
order of the soliton lattice destroys. By overcoming the 
nucleation energy barrier of the order parameter vortices, 
a first order phase transition from easy plane spiral state 
to a saturated ferromagnetic state occurs at h = h c %, 
depends on J s , and Q. This transition turns out to a 
continuous cross-over if the initial state of the system at 
h = h c i holds the zero spin-supercurrent (a uniform easy 
plane with Q = 0). The possibility of the existance of 
the superfluid phase and crossing onto a saturated ferro- 
magnetic state (the dependence of h c i on Q, and J s ) are 
discussed in some detail. 

The zero-temperature quantum ferrimagnetic chain is 
defined by H = JJ2<ij> ' s i ~ + s «) where 

S =/= s. The low energy physics of the quantum ferri- 
magnets in the presence of the external magnetic field 
can be obtained by a 1-dimensional NLctM. The external 
magnetic field couples with n, the Neel order parameter 
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C h = iAfoA(n) • d T n + C NLaM +h + Aop - M h ■ n, (1) 
where 

£riL CT M+h = 7r (v s (d x n) 2 + — [2ih + n x d T n} 2 ] . (2) 

The first term in Eq. ([!]) is the usual (dynamical) Berry's 
phase of a quantum ferromagnet, with Mo = \S — s|/&o 
the magnetization per unit cell (pair of sites). Here 
£top is the topological term, g = 4/(s + S), and v s — 
4a JsS/(s + S) ||. It is this term that results in the 
ferromagnetic branch of the spin waves and corresponds 
to the trajectory of spin over a closed orbit on the unit 
sphere in the presence of a unit magnetic monopole at the 
center. The contribution of the first term in Ch is equiv- 
alent to the area enclosed by this trajectory and since ei- 
ther of the smaller or the larger enclosed areas on the unit 
sphere must lead to the same Berry's phase, the magnetic 
moment per unit cell, i.e. 2Aioao, must be quantized with 
integral values |?|]. Setting Mo — in action (|l|) follows 
to the usual 0(3) NLcrM in the magnetic field ||], appli- 
cable to the Heisenberg antiferromagnets. Similar to zero 
mag netic field g we can find the spin-wave modes. At 
h = the ferrimagnetic spin- waves consist of both gapless 
(ferromagnetic) and gapped (antiferromagnetic) modes, 
and at T = the low (high) energy physics of quantum 
ferrimagnets is effectively like that of a ferromagnet (anti- 
ferromagnet) which is formed by the chain of (dimerized) 
unit cells with magnetic moment Mq(= \S — s\). Apply- 
ing an external magnetic field, h, develops a gapped ferro- 
magnetic spin- waves; in which case the energy cost for the 
ferromagnetic transitions is proportional to the Zeeman 
splitting factor. Unlike to ferromagnetic mode, the effect 
of the external magnetic field is to suppress the antifer- 
romagnetic gap. Clearly, this reveals similarities between 
the ferrimagnets and the integer spin Heisenberg antifer- 
romagnets in magnetic field [|9| . The ground state of the 
fcrrimagnet corresponds to the staggered configuration of 
spins, unless h > h c \ = 2J\S — s\. At this point the stag- 
gered state becomes unstable against the non-collinear 
spin-flop phase (the partially polarized state) of the spins 
when the spectrum becomes soft at k = 0. When the ex- 
ternal magnetic field exceeds h C 2, the system will be in a 
saturated ferromagnetic phase, with a quantized magne- 
tization per unit cell. The h C 2 = 2 J(S + s) is obtained 
by using the dispersion relation of the spin waves based 
on the fully polarized state of the ferrimagnets. It is the 
lower-bound of the external magnetic field, in the sense 
that the spin waves of the fully polarized state become 
soft at k = 0. 

The outline of this paper is as follows : In next section 
we will present the numerical computations of Lanczos 
method to compare the behavior of the correlation func- 
tions on the plateau and between two plateaux. We also 
derive the critical exponents, and the effective Hamilto- 
nian for the latter region by using a Quantum Renormal- 



ization Group (QRG) approach. The main part of this 
paper, the spiral (superfluid) state of the quantum fer- 
rimagnets is presented in section III. The details of the 
collective modes in the intermediate magnetic phase is 
given as the appendix in Sec. IV. 

II. NUMERICAL RESULTS 

In order to have more accurate physical picture, we 
present the numerical results, by applying a Lanczos 
method to each 5 2 -sector of Hilbert space from S z = 
f (S-s) to S, = f (S+s). Here N(= 20) is the number 
of sites that is used in the exact diagonalization. 

A curve for magnetization vs. magnetic field for a 
chain of (1/2, l)-ferrimagnet has been given in Ref. [fiol . 
An extrapolation to N — > oo on exact diagonalization 
calculation shows there are two plateaux of magnetiza- 
tion below h c i — Ao — 1.7589J and above h C 2 = 3J, 
associated with the magnetization m = = 1/3 and 
m = 1. Here A is the energy gap between the Ferro- 
and Antiferromagnetic spin wave modes of this model. 
In Fig. [j] the numerical results of the in-plane spin-spin 
correlations is presented for a point on the plateau (at 
m = 1/3) and some intermediate points between two 
plateaux, i.e., m = 2/5,8/15,2/3 and 4/5. This exhibits 
the in-plane spin-spin correlation functions within the in- 
termediate magnetic field region which falls off as power 
law, and manifests the critical behavior. 
Thus the in-plane correlation is expected to have the 
asymptotic form (S x (0)S x (r)) ~ r~ v . For instance we 
have calculated this exponent for m=2/3 by using data 
of exact diagonalization of a chain with length N = 20. 
Since there are correlations between different types of 
sublattices we find an exponent for each case. For cor- 
relations on sublattice A (Fig.|l|-(a)) and m=2/3 we ob- 
tained r\ = 0.44 ± 0.01 and for sublattice B (Fig.|-(b)) 
we have r\ = 0.42 ± 0.04. A similar behavior is seen for 
the correlations of different sublattices A and B (Fig.|l|- 
(c)) where r\ = 0.47 ± 0.01. To have more qualitative 
picture of the transient region between two plateaux, we 
can use a quantum renormalization group. To perform 
this, we choose two adjacent spins (S = l,s = 1/2) 
as the building block. The block Hamiltonian is then 
Hb = JS-s — h(S z + s z ). For low field limit (e.g. h < h c {) 
the lowest lying states of Hb are a spin- 1/2 doublet. This 
comes out with an effective Hamiltonian of ferromagnetic 
Heisenberg chain with S = 1/2 in magnetic field, corre- 
sponds to the m = 1/3 plateau. For high field limit 
h > 3J/2 where two states \S T = 1/2, Sj. = 1/2), and 
\St = 3/2,51, = 3/2) are nearly degenerate, we arrive 
with a spin-1/2 antiferromagnetic XX Z + h Hamilto- 
nian. The effective Hamiltonian is 

N/2 N/2 

tfcff = "J E( T « ' T n + 1 + A« +1 ) ~h'J2< (3) 
n— 1 n—1 
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FIG. 1. The log- log plot of in- plane spin-spin correlation 
function is shown within the plateau (m = 1/3) and within 
the partially spin polarized state (m = 2/5,8/15,2/3,4/5). 
The former shows the exponential decay while in the latter 
region the power law behavior of correlation function is seen. 
The correlation on sublattice A (spin-1, spin-1) is shown in 
(a), on sublattice B (spin-1/2, spinl/2) in (b) and between 
the two sublattices A (spin-1) and B (spin- 1/2) in (c). In (c) 
spin-1 is fixed and the position of spin-1/2 is running to show 
the correlations. In each part an inset represents the corre- 
lation for m = 1/3 and 2/3 in a normal scale. The former 
approaches to zero exponentially, and the latter falls off alge- 
braically. 

where r is spin-1/2 operator, A = 1/3, h' = |/i — ^J. 
It is known such model Hamiltonian in the presence of 
a magnetic field (h) has a critical line which separates 
the partially polarized phase from the fully polarized one 
p| . The magnetization of the ferrimagnetic chain (to) 
is related to the magnetization of XX Z + h model by 
m = 1 + mxxz+h which leads to the same behavior 
as the numerical results but with some discrepancies for 
the critical fields (h c i = 1.5J , h C 2 — 3 J). However 
one may continue the quantum RG procedure for the 
XXZ + h model and obtain that the RG flow for the par- 



tially polarized phase (between the two plateaux) goes to 
the isotropic XY fixed point J12|. 



III. SUPERFLUID PHASE 

By generalizing the NLerM+h we propose a phe- 
nomenological field theory which can present the long 
wave length limit of the quantum ferrimagnets within 
h > h c i 



E = dx 
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:pf(d x n ± Y + -pl(d x n z y + (3ni 



(4) 



where we applied the spin coherent state formalism in 
Eq.(|j). Here h* = h — h r \ is the effective magnetic field. 
pjr and p z s are in-plane and out of plane spin stiffness. 
(3 > gives an easy plane anisotropy. To study a su- 
pcrfluid phase, it is convenient to introduce the following 
variational solutions n(x) = (sin 9 cos tp, sin 9 sin tp, cos 9) , 
where 9 = 9(x) and tp — Qx. Here Q is a spiral 
state wave vector responsible for the superfluid phase 
slip, and a = 2tt/Q is the soliton lattice spacing, i.e., 
the phase tp changes by 2tt along the chain. Clearly 
h z = cos 9 is the uniform solution (6 is a constant) and 
E = pjrQ 2 (l - h 2 z )/2 + j3h 2 z - h*n z is the correspond- 
ing energy per unit length. Note n z is the momentum 
density conjugate to field variable tp, and the Hamilto- 
nian (||) gives a linearly dispersing collective mode, asso- 
ciated with the U(l) symmetry breaking phase, i.e., the 
superfluid phase. It is easy to check how the classical 
solutions and the fluctuating modes can be derived by 
Eq.(^). For example dE/dh z = leads to n z — h* /K zz 
where K zz = d 2 E/dn 2 z — 2(3 — pjrQ 2 is the energy gap 
of the out-of-plane modes at zero wave vector. Then 
J s = (l/h)dE/dQ — n\pjrQ/Ti is the gauge invariance 
topological spin-supercurrent density carried along the 
x-direction. Furthermore, it is necessary to study the 
stability of the superfluid phase against the quantum 
fluctuations. This governs with the in-plane fluctuations 
K vv {k) = 2p-jrn 2 L k 2 , and the out of plane fluctuations 



K, 



,(k)=pfni(k 2 + Q 2 )-pfhiQ 2 

+p z s fi\k 2 + 2(3{n\ - h 2 z ) + h*h z . 



(5) 



As one can see K z 



K zz {k = 0) at h* = 0. The 



gapless linear superfluid mode can be obtained by Tito — 
2^/KtpspKzz. The details of the collective modes calcula- 
tion is presented in Appendix. Through out this formula- 
tion, the zero temperature phase diagram of this system 
can be obtained. At h* — where h z — the out of plane 
instability occurs at Q c = At this point the en- 

ergy gap of the zz-modes vanishes and the 0(3) symme- 
try of the Hamiltonian is restored (the dispersion relation 
of the collective modes become imaginary at small k) . It 
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follows the transition to a uniform solution must occurs 
at this point since iri(S 2 ) — where the linear solution 
can be considered as a map from the compactified phys- 
ical space to the equator of the order parameter space 
S 2 . Such instability in order parameter space is induced 
by the local fluctuations of the order parameter, out of 
the equator toward the north pole of the order parameter 
space S 2 . When J s is close to its critical values (J sc ), the 
local fluctuations are very strong, and it is likely the order 
parameter pass the north pole, and the phase p becomes 
singular, i.e., a vortex in p is nucleated, and the phase 
winding is lost. The effect of h* is pushing the spins to be 
aligned along z-axis, e.g., n z = 1 at large enough h*. We 
notice if h > h c i, system undergoes onto the saturated 
ferromagnetic phase by increasing the spin-supercurrent 
density. This happens at Q* = \f{2j3 - h*)/pj- < Q c , 
corresponding to J s = 0, before the outset of the easy- 
plane uniform state. To study the nucleating mechanism 
of the vortices, we implement a mechanical analogy to 
the classical field theory. This formalism has been de- 
veloped by Langer, and Ambegaokar [|| for one dimen- 
sional superconductors. To avoid the complexity let us 
neglect p z s in energy functional (|i|) from now on. To 
start this calculation, we make the following transforma- 
tions n x {x) = f (x) cos <p(x) , n y (x) = f (x) sin ip(x), and 
n z {x) = y/l — f 2 (x) (where / = n±), and 



E[f] 



{d x ff + j 2 



+P(i-f 2 )-h*sr~p) 



(6) 



where L is the momentum conjugate to <p and it is pro- 
portional to the spin-supercurrent density J s , because 
d x p = L/f 2 is the solution of the 5E[f,ip]/5(p = 0. It is 
straightforward to show how SE[f]/Sf = yields 



x — Xq 



/(*) 



df 



f(x ) 2 v / E cS - Uctiif) ' 



f{x) dx 



ip{x) - tp(x ) = L I — 
f(x ) J 



(7a) 



(7b) 



Ues(f) 



L 



2/2 



+ [3(l-f 2 )-h*^/T~f 2 ~. (7c) 



Uniform solutions (f — f — etc) arc one of the solu- 
tions of Eqs. (0). It leads to p — Qx and L = Qf 2 . 
The energy potential associated with the uniform solu- 
tions are depicted in Fig. || where — U e g(f) vs. / is 
plotted for h* = J for different Q's, or equivalently for 
different spin-supercurrent J s . Here pjr = 2J/3, and 
(3 = J. The effective potential, —U e s, has one minimum 
at /= 1 - [h*/(2/3 - P iQ 2 )] 2 ± if Q < Q*. 
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FIG. 2. — U e ff(f), the effective potential of the spiral state 
vs. / is plotted for h* = J. The minimum energy solution at 
Q — Q* is / = where the spiral state destroys. 



This minimum energy solution is disappeared if Q > 
Q*(= \J (2/3 — h*)p s ), and is replaced by another min- 
imum at / = 0. At Q* a first order phase tran- 
sition from / 7^ onto / = takes place by nu- 
cleating of a vortex. One can follow the transition 
between these minimum energy solutions by changing 
the parameters, Q and h* , as illustrated in Fig. |^. 
The saddle point solutions can be derived by setting 
d 2 U c s/df 2 — at the extremum of U e g(f). It follows 
at Q* = y/{2/3 - h*)/pj- (or equivalently J s = 0), the 
local minima at finite / is disappeared and n z = I be- 
comes a unique minima of U e g. Moreover, the cross- 
over from the uniform state [a broken U(l) symmetry 
state with zero spin-supercurrent] to a saturated ferro- 
magnetic state occurs at h* — 2(3 — 2 J (Q* = 0), con- 
sistent with the exact diagonalization results (h C 2 = 3J). 
Obviously at h = h c \ and Q = Q c (= y/2(3/p s ), the 
sign of curvature of —U e g changes, and therefore Q c can 
be considered as the saddle point of —U c g, at which the 
spiral state becomes unstable against the hard-axis fluc- 
tuations and a transition to a uniform solution occurs. 
This is seen if the spin-supercurrent J s arrives to its crit- 
ical value, J sc . As was mentioned earlier, the mechanism 
for such transition is nucleating of a vortex in the order 
parameter space. For a given Q and h*, it is straight- 
forward to show E(Q) = p^Q 2 /2 - h* 2 /{2(3 - p^Q 2 ). 
Recalling J s = (1 / %)dE{Q) / 'dQ leads to J s = pjf 2 Q/h 
(J s = p-jrL/K) which is the spin-supercurrent density and 
n z = h*/(2f3 — pjrQ 2 ) is the classical solution of Eq.(||). 
The dependence of L(= Q[l - h* 2 /{2(3 - PsQ 2 ) 2 }) with 
respect to Q for various h* is illustrated in Fig. [| It 
follows that J s vanishes at Q* = \J {2(3 — h*)/pj- where 
the system crosses to the saturated ferromagnetic phase. 
The maximum spin-supercurrent J* (< J sc ), which can 
pass through the system at finite wave vector Q (and 
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h* ^ 0), is a monotonically decreasing function of h* . 



dtpjq) 
dt 



6H f 



5^n z (-q) 



= T K zz(q)n z (q), 



(9b) 




FIG. 3. L vs. Q is plotted for different h* . For any h* there 
is a maximum spin-supercurrent J s = J* (< J sc ) which can 
pass through the system at wave vector Q' . J| is a mono- 
tonically decreasing function of h* . At Q — Q* 7^ 0, the 
spin-supercurrent density vanishes [L = 0). 

In conclusion, we have predicted the intermediate mag- 
netic field region of fcrrimagnets can support the dissi- 
pationlcss flow of the spin-supercurrents, which may be 
observed by the advance techniques of spintronics. The 
idea of the collective spin-supercurrent transport which 
has been presented in this paper is completely general, 
and can be applicable to any 1-dimensional spin system 
with easy plane anisotropy. The intermediate magnetic 
field phase of the ferrimagnets is one example. Other ex- 
ample of this kind is the isospm-supercurrent transport 
of the quantum Hall bars in the bilayer electron systems. 



where n z is the momentum density associated with the 
field variable ip. In the spin coherent representation, n z 
and n+ are given by —ihd/dip, and e i<p respectively. Then 
it is easy to check the transformation n z — > dtip can be 
obtained by canonical quantization, i.e., n z can be con- 
sidered as the momentum density conjugate to field vari- 
able ip. Making the derivative with respect to time we 
find the equation of motion: fi z = — {2/K)K tptp {q)ip(q) = 
— {2/K) 2 K tpv {q)K zz {q)n z , where ip = dip/dt. This clearly 
gives 



huj = 2 J K v Jq)K zz (q) 



(10) 



Eqs.(|9a-9b) are similar to the coupled Josephson junction 
relations in superconductivity. 
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IV. APPENDIX [1] 

Here we present the detailed calculation of the collec- 
tive modes in the intermediate magnetic field phase. The PI 
effective Hamiltonian in terms of the in-plane and the out 
of plane fluctuations is given by 

H f [ip,n z ] = ^^2<p(-q)K 9v (q)ip{q) 

+ \^Tn z {-q)K zz (q)n z {q) : (8) 



where K w = S 2 E/dip 2 and K zz = S 2 E/Sn 2 z . The equa- 
tion of motion for the field variables ip, and n z can be 
obtained by the Hamilton equations (p = —5H/5q and 
q = 6H/6p) 

dn z (q) _ SH f 2 



dt 



S'M-q) 



— K 



(9)^(9): 
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